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LINEARIZED  COMPRESSIBLE-FLOW  THEORY 

FOR  SONIC  FLIGHT  SPEEDS 

By  Max.  A.  Heaslet,  Harvard  Lomax,  and 
John  R.  Spreiter 


SUMMARY 

/ 

The  partial  differential  equation  for  the  perturbation  velocity 
potential  is  examined  for  free-stream  Mach  numbers  close  to  and  equal 
to  one.  It  is  found  that,  under  the  assumptions  of  linearized  theory, 
solutions  can  be  found  consistent  with  the  theory  for  lifting-surface 
problems  both  in  stationary  three-dimensional  flow  and  in  unsteady 
two-dimensional  flow.  Several  examples  are  solved  including  a  three- 
dimensional  swept-back  wing  and  a  two-dimensional  harmonically 
oscillating  wing,  both  for  a  free-stream  Mach  number  equal  to  one. 


INTRODUCTION 

Much  of  the  recent  progress  in  the  theoretical  analysis  of 
compressible— flow  fields  is  attributable  to  the  successful  applica¬ 
tion  of  linearization  methods.  Although  the  basic  assumptions  used 
In  conventional  linearized  theory  appear  at  first  glance  to  be  highly 
restrictive,  it  has  been  found  that,  just  as  in  the  analogous  case  of 
thin— airfoil  theory  for  incompressible  flow,  the  methods  have  many 
fields  of  utilization  adequate  for  most  engineering  purposes.  Since 
the  basic  methods  are  so  well  known  and  depend  on  such  relatively 
simple  mathematical  tools,  it  appears  obvious  that  the  range  of 
applicability  of  the  theory  should  be  explored  completely*,  Such  is 
the  purpose  of  the  present  report.  It  has  been  more  or  less  tacitly 
presumed  in  the  past  that  such  applications  cannot  treat  cases  for 
which  the  flight  velocity  is  near  the  speed  of  sound.  In  the  study 
of  two-dimensional  steady— state  problems  in  airfoil  theory,  this 
presumption  is  certainly  true.  The  Prandtl-Glauert  and  Ackeret  rules 
for  variation  of  pressure  coefficient  with  free-stream  Mach  number 
in  the  subsonic  and  supersonic  regimes,  respectively,  are  clearly 
invalid  for  Mach  numbers  near  one,  since  perturbation  velocities 
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become  arbitrarily  large.  In  this  case ,  linearized  theory  therefore 
predicts  its  inability  to  treat  such  problems.  On  the  other  hand, 
if  linear  methods  are  applied  to  nonstationary  two-dimensional 
airfoil  and  particular  steady— state,  three-dimensional,  lifting- 
surface  problems  at  sonic  speeds,  a  consistent  theory  results  since 
solutions  are  found  which  yield  perturbation  velocities  of  the  same 
order  of  magnitude  as  those  calculated  for  free-stream  Mach  numbers 
of,  say,  0.6  or  1.5. 

Unfortunately,  arbitrary  thickness  distributions  at  sonic  speeds 
cannot  be  studied  by  linear  theory  in  the  steady  state  since,  in 
general,  the  theory  predicts  infinite  pressure  differences  between 
the  wing  surface  and  infinity.  In  the  particular  case  of  a  yawed, 
symmetrical  wing  of  infinite  aspect  ratio,  the  results  are,  however, 
again  consistent  with  the  theory  and  yield  pressure  distributions 
which  are  the  same  as  those  determined  by  using  only  the  component 
of  free— stream  velocity  normal  to  the  leading  edge.  The  derivation 
of  this  latter  result  for  a  free— stream  Mach  number  of  one  will  be 
given. 

The  difficulty  of  not  being  able  to  include  thickness  effects 
in  general,  together  with  the  uncertainty  of  the  magnitude  of  the 
viscous  effects,  leaves  the  question  as  to  the  limitations  of  such 
a  linear  theory  in  application  to  practical  wing  shapes.  Such  a 
question  can  certainly  not  be  resolved  by  mathematical  reasoning 
alone.  The  extent  to  which  the  fluid  medium  can  be  idealized  at 
these  speeds  is  left,  for  the  time  being,  unsettled  and  it  remains 
for  experiment  to  determine  whether  the  consistent  mathematical 
results  which  are  obtained  from  the  linearized  equations  provide 
reasonably  exact  predictions.  In  this  connection,  it  should  be 
mentioned  that  the  few  experimental  results  available  for  the  total 
'  lift  on  thin  triangular  wings  at  Mach  numbers  near  one  tend  to 

confirm  the  theory.  But  even  if  more  detailed  experimental  results 
indicate  that  further  refinements  are  necessary,  there  is  still 
little  doubt  but  that  the  linear  potential  solutions  will  provide  a 
valuable  basis  for  more  exact  extensions  of  theory. 

The  present  report  is  divided  into  three  parts.  In  the  first 
part,  the  linearization  of  the  partial  differential  equation  for 
the  velocity  potential  is  carried  out  in  some  detail  for  steady- 
state  conditions.  A  by-product  of  this  derivation  is  the  nonlinear 
form  of  the  equation  for  two-dimensional  flow  which  was  used  by 
von  Kaxmaxi  (reference  l)  to  determine  his  similarity  rules  for  tran¬ 
sonic  flow.  The  equation  for  unsteady  two-dimensional  flow  based  on 
the  same  assumptions  is  also  given.  The  second  part  of  the  report 
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is  restricted  to  two-dimensional  unsteady  problems  for  values  of 
Mach  number  near  one.  The  principal  contribution  of  this  section  is 
the  evaluation  of  the  change  with  time  of  the  pressure  distribution 
over  an  airfoil  starting  suddenly  from  rest  at  a  speed  close  to  that 
of  sound.  Such  an  idealized  problem  involves  a  step  function  in 
velocity  in  which  the  airfoil  has  zero  velocity  for  all  negative 
values  and  near  sonic  velocity  for  all  positive  values  of  time.  From 
these  results  the  initial  build-up  of  lift  can  be  calculated  for  Mach 
numbers  near  one,  although  the  eventual  value  of  the  lift  cannot  be 
found  by  linear  methods.  Further  application  can  also  be  made  to 
problems  in  flutter  and  gust  loads.  The  third  part  of  the  report 
treats  the  steady— state  three-dimensional  problem.  Both  lifting 
surfaces  and  symmetrical  nonlifting  wings  are  considered  and  it  is 
seen  that  in  the  former  case  consistent  solutions  are  obtained  by 
particularly  simple  means.  These  solutions  represent  the  limiting 
case  of  both  subsonic  and  supersonic  lifting-surface  theory  and 
give,  for  example,  the  same  value  of  lift— curve  slope  at  the  speed 
of  sound  that  was  obtained  for  the  supersonic  triangular  wing  by 
Stewart  (reference  2). 

A  list  of  symbols  is  given  in  the  appendix. 


PART  I  -  THE  LINEARIZED  EQUATIONS  OF  MOTION 
Steady  State 

The  nonlinear  partial  differential  equation  satisfied  by  the 
velocity  potential  $  of  an  isentropic  flow  field  can  be  expressed 
in  the  form 


$xx 


—  ®yz  £y$z - —  ®zx  ^z^x  — —  ®xy  ®x®y  =  0 


(1) 


where  the  subscript  notation  is  used  to  indicate  differentiation  and 
a  is  the  local  speed  of  sound  given  by  the  relation 


In  this  latter  equation  Vq  and  Mq  are,  respectively,  velocity 
and  Mach  number  of  the  free  stream,  y  is  the  ratio  of  specific 
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heats  (for  air,  7=1.4),  and  V  is  local  velocity. 

Introducing  the  perturbation  velocity  potential  cp,  where 

cp  =  -Vo  x  +  $  (3) 

it  is  possible  to  express  equation  (l)  in  terms  of  the  derivatives 
of  cp  and  the  parameters  Mq  and  V0.  To  begin  the  linearization 
of  the  resulting  equation,  the  coefficients  of  the  second  ordered 
derivatives  of  $  are  expanded  in  Maclaurin  series  with  ascending 


powers  of 


The  convergence  is  assured  provided 


u2+v2+w2  ^ 


<1 


(4) 


or,  in  a  slightly  modified  form,  provided 


Pa  2 
7-1 


If  the  assumption  is  now  made  that  -p-,  -jf->  «  1  so  that 

vo  vo  vo 

second  and  higher  powers  in  the  perturbation  velocities  can  be 
neglected  in  comparison  with  one,  the  partial  differential  equation 
can  be  simplified  to  the  form 


* 


* 


1  —  Mq2 


1  +  —  +  (7-1)  M  2 


u 

V„ 


j-+  ^yy  +  4*zz 


—  2  <T>xz  —  Mq2  —  2  'Pxy  —  Mq2  =  0  (6) 

V  V 

vo  o 

From  this  equation  all  the  succeeding  expressions  will  be  derived. 

Two-  and  three-dimensional  linear  equations,  Mq  £  1.  -  Since 

equation  (6)  is  obviously  nonlinear,  additional  assumptions  must  be 
made  to  reduce  it  to  a  linear  form.  Clearly,  these  assumptions  must 
involve  the  relative  magnitudes  of  all  the  terms  in  order  to  determine 
which  ones  may  be  neglected.  Perhaps  one  of  the  least  restrictive 

set  of  conditions  is  that:  * 
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(a)  The  ratios  of  the  perturbation  velocities  to 
the  free— stream  velocity  are  small  enough  to  be 
neglected  when  compared  to  one. 

(b)  The  velocity  gradients  at  a  given  point  of 
the  flow  field  are  all  of  similar  magnitude. 

With  the  aid  of  these  assumptions,  it  follows  that,  to  the 
order  of  the  approximations  made,  the  perturbation  velocity  potential 
cp  satifies  the  well-known  linear  equation 

(1  -M02)  <PXX  +?yy  +<P2Z  =  0  (7) 

In  the  case  of  two-dimensional  flow,  the  equation  is  independent  of 
y  and  thus  may  be  written  in  the  form 

(1  -Mo*>  *0  +tfW  0  (8) 


Two—  and  three-dimensional  nonlinear  equations,  Mn  =  1.  — 

The  study  of  equation  (8)  in  both  subsonic  and  supersonic  flow  has 
shown  that  for  arbitrary  lifting  surfaces  or  symmetrical  nonlifting 
airfoils  the  value  of  the  induced  velocity  u  on  the  surface  of  a 

fixed  geometric  configuration  is  proportional  to  ( |  1  -  mo2 | )~X^ . 

In  all  airfoil  problems,  the  value  of  u  becomes  infinitely  large 
as  Mq  approaches  one,  either  from  above  or  below,  and  the  basic 
assumptions  are  thus  violated.  Such  a  difficulty  led  Oswatitsch 
and  Wieghardt  (reference  3)  and  Sauer  (reference  4)  to  abandon  the 
restriction  of  linearity  and  to  seek  a  more  exact  equation  at  M0  =  1. 
Retaining  the  assumptions  underlying  equation  (6)  and  setting  ° 

V0  =  a*  where  a*  is  the  critical  speed  of  sound,  it  follows  that 
at  M0  =  1  the  perturbation  velocity  potential  satisfies  the  equation 


(7+1)  2 

a*  9X  ‘Pxx  (Pzz+  a*  9z'cPXz 


=  0 


Since  cpx  is  much  larger  than  <p2  as  the  Mach  number  approaches 
one,  equation  (9)  may  be  further  simplified  to 


(7+1) 

a* 


<P, 


<P. 


xx 


cp 


zz 


=  0 


(10) 


If,  in  three  dimensions,  the  perturbation  velocities  do  not 
remain  small,  equation  (6)  again  supplies  the  necessary  form  of  the 
differential  equation  at  MQ  =  1.  From  the  relation  Y  =  a*,  the 
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required  expression  is 


(7+1) 

a* 


^xx 


y 


(n) 


Two—  and  three-dimensional  linear  equations,  MQ  =  1*  — 

Equation  (10)  has  been  used  by  von  Karmaa  (reference  l)  to  establish 
similarity  rules  for  two-dimensional  transonic  flow  and  is  the  basis 
for  work  continuing  at  the  present  time.  (See  also  reference  5») 

If  at  MQ  «  1  the  assumptions  made  in  the  linearization  process  still 
hold,  it  follows  from  either  equation  (10)  or  (8)  that  the  differential 
equation  reduces  to  the  form 

<pzz  =  0  (12) 


It  is  possible,  however,  to  predict  independently  from  this  relation 
that  linearized  methods  cannot  be  applied  to  the  calculation  of 
arbitrary  airfoil  pressure  distributions.  The  range  of  applicability 
of  such  an  equation  is  thus  almost  nonexistent.  On  the  other  hand, 
the  linearized  form  of  equation  (ll)  or  (7)  at  M0  =  1  is 

<Pyy  +  Vzz  =  0  (13) 

and  from  this  equation  a  class  of  nontrivial  solutions  can  be 
obtained  for  particular  boundary  conditions.  Both  equations  are 
of  parabolic  form  in  the  number  of  dimensions  for  which  they  are 
defined.  In  the  present  report,  formal  solutions  satisfying  the 
imposed  conditions  will  be  obtained  in  three  dimensions  for  flat 
lifting  surfaces  with  swept-back  leading  edges  and  for  an  infinitely 
long,  symmetrical,  swept— back  wing. 


V 


Unsteady  State 


The  derivation  of  the  steady-state  equations  for  the  velocity 
potential  was  developed  in  some  detail  because  of  the  various  results 
to  be  obtained.  Similar  methods  can  be  used  when  unsteady  conditions 
are  to  be  considered,  the  differential  equation  for  the  velocity 
potential  being  now  in  the  form 


~  ^  ($t*t»  +  2$x  ®xt*  +  2<*y  $ytf  +  2$z  ®zt *) 

-  ®yz  ^z  -  ®zx  *z  $x  -  §f  ®xy  ®x  ®y  =  0 


4 


(14) 


where  tf  represents  time.  The  details  of  the  derivation  can, 
however,  be  avoided  by  referring  directly  to  the  equation  satisfied 
by  the  velocity  potential  for  the  propagation  of  sound  waves  of 
small  amplitude.  (See  reference  6,  p.  i-92.)  In  this  form  of  the 
equation  the  Cartesian  coordinate  system  x,  y,  z  is  assumed  fixed 
in  the  medium  so  that  free— stream  velocity  is  zero,  while  the  wing, 
which  moves  in  the  direction  of  the  negative  x  axis  with  velocity 
VQ,  generates  small  pressure  disturbances.  As  a  consequence,  the 
velocity  potential  of  the  field  satisfies  the  well-known  wave  equation 
in  three  space  dimensions: 

—  “  CPxx  “  Tyy  “  ^22  =  0  (15) 

Equation  (15)  is  reducible  to  canonical  form  by  means  of  the 
relation 


t  =  a0tf 

and  the  three-dimensional  form  of  the  equation  is  therefore 

<Ptt  ~<Pzx  _cPyy  _cPzz  =  0  (l6) 

while  in  the  two-dimensional  case  independence  with  respect  to  y 
yields 

-  V  -  V  * 0  <17> 


PART  II.  TWO-DIMENSIONAL  LINEAR  PROBLEMS  FOR  Mo  NEAR  ONE 
Unsteady  State,  Mq  =  1 

It  was  pointed  out  in  the  derivation  of  equation  (12)  that  the 
linear  equation  for  the  velocity  potential  is  not  applicable  to  air¬ 
foil  problems  in  either  the  subsonic  or  supersonic  regimes  for  Mq 
near  one.  The  possibility  still  remains,  however,  of  analyzing 
unsteady  flows  during  the  period  in  which  the  perturbation  velocities 
remain  small.  As  an  example  of  such  a  problem,  consider  the  case  of 
a  flat  lifting  surface  at  a  small  angle  of  attack  a  starting  from 
rest  at  a  velocity  V0  near  the  speed  of  sound.  The  perturbation 
potential  for  such  a  motion  is  equivalent  to  the  change  in  potential 
brought  about  by  an  abrupt  change  a  in  angle  of  attack  of  an  air¬ 
foil  flying  in  a  steady— state  condition  at  velocity  equal  to  V0„ 
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The  change  in  load  distribution  brought  about  by  this  maneuver  is 
the  so-called  indicial  load  distribution  and  a  knowledge  of  such 
indicial  functions  is  important  fdr  applications  using  operational 
calculus. 

- characteristic  traces  Figures  1(a)  and  1(b) 

furnish  an  insight  into 
X  the  nature  of  the  boundary 

conditions  for  airfoils 
traveling,  respectively, 
at  supersonic  and  subsonic 
speeds.  The  chord  is  ini¬ 
tially  on  the  x  axis  with 
leading  edge  at  the  origin 
and  trailing  edge  at 

(a)  Supersonic  wing.  x  =  c0.  With  increasing 

time,  the  wing  section 
X  travels  in  the  negative 

x  direction  and  sweeps 
out  a  portion  of  the  xt 
plane  as  shown  in  the 
figures  (indicated  by 
shaded  areas).  Through¬ 
out  this  part  of  the 
plane  the  boundary  condi— 

f  .  t  tions  require  that  the 

(b)  Subsonic  wing.  induced  vertical  velocity 

is  — V0a,  while  elsewhere 
Figure  1.-  Boundary  conditions  for  on  the  plane  the  lnduced 

two-dimensional  unsteady-lift  velocities  are  continuous 

problem.  functions  of  z.  In  fig¬ 

ure  2,  sketches  of  the 

airfoil  in  the  supersonic  and  subsonic  cases  are  shown  together 
with  indications  of  the  manner  in  which  the  disturbance  field 
spreads.  These  wings  are  presented  in  xyz  space  with  time  as 
a  parameter  so  that  their  coordinate  system  is  not  to  be  confused 
with  that  of  figure  1.  The  airfoils  are  traveling  from  right  to 
left  at  Mach  numbers  of  0.8  in  the  subsonic  case  and  1.2  in  the 
supersonic  case,  and  for  a  time  corresponding  to  that  required 
for  the  wing  to  travel  a  distance  of  one— third  chord  length.  At 
t  =  0  cylindrical  waves  are  induced  at  each  disturbance  point,  that 
is,  at  each  point  of  the  chord.  These  waves  expand  radially  at  the 
center  of  the  expanding  waves  moves  relative  to  the  initial  disturb¬ 
ance  point.  At  a  given  instant  in  time  the  entire  disturbance  region 
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of  the  airfoil  is  contained  within  the  closed  surfaces  shown  in  the 
figures,  the  outer  surfaces  corresponding  to  the  largest  values  of 


Figure  2.-  Sketch  showing  extent  of  disturbance  fields  after 
travel  of  one— third  chord  length. 

time.  In  the  supersonic  case,  the  pressure  distribution  over  the 
wing  reaches  a  steady-state  value  as  soon  as  the  wing  moves  ahead  of 
the  expanding  cylindrical  wave  produced  at  t  =  0  by  the  leading 
edge.  In  the  subsonic  case,  the  wing  never  leaves  the  disturbance 
field  of  the  cylinders  and,  as  will  be  seen  later,  the  steady— state 
pressure  distribution  is  approached  asymptotically. 

It  is  apparent  from  equation  (17)  that  the  characteristic  cones 
have  semivertex  angles  equal  to  45°  and  that  the  cones  with  vertices 
on  the  xt  plane  have  traces  with  slopes  equal  to  ±1.  These  cones 
determine  the  upstream  boundary  of  the  field  of  influence  of  the 
vertex  point  and  their  cross  sections  in  the  plane  t  =  constant 
are  the  disturbance  regions  of  the  cylindrical  waves  arising  at  the 
vertex.  Thus,  perturbations  in  pressure  produced  initially  at  the 
leading  edge  of  the  wing  section  are  confined  at  later  time,  in  the 
modified  coordinate  system,  to  the  cone  with  vertex  at  the  origin 
and  traces  x  =  ±t. 
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The  solution  of  equation  (17)  for  boundary  values  of  the  type 
under  discussion  has  been  indicated  in  reference  7  through  considera¬ 
tion  of  an  analogue  problem  in  supersonic  lifting— surface  theory. 
Thus,  the  shaded  areas  in  figure  s;  l(a)  and  1(b)  are  thought  of  as 
swept— forward  lifting  surfaces  situated  in  a  stream  directed  along 
the  positive  t  axis  at  a  Mach  number  M o  =  ,y/2.  The  boundary  values 
remain  the  same;  that  is,  <Pz  =  w  =  — V0a  o n  the  wing  and  ‘Py*  ^z 

are  continuous  functions  of  z  elsewhere  in  the  xt  plane.  In 
lifting— surface  terminology,  the  unsteady  case  for  supersonic  speed 
becomes  a  wing  with  supersonic  leading  edge,  while  the  case  indicated 
in  figure  l(b)  involves  a  subsonic  leading  edge. 

The  solution  for  the  wing  traveling  at  supersonic  speed  has  been 
given  in  reference  7  in  a  form  valid  for  all  Mach  numbers  greater 
than  or  equal  to  one.  The  expressions  for  load  coefficient  Ap/q, 
where 

4P  _  Pz^u 

<1  ~  In  V2 

2P°V° 


differ  analytically  in  various  regions  of  the  xt  plane.  These 
expressions  are: 

Region  A  (between  lines  x  =  -M0t,  x  =  -t) 


$ 


(l8a) 


Region  B  (between  lines  x  =-t,  x  =  t,  and  x  =  c0-MQt) 


£e  = 
q 


4a, 


—  arc  cos 

JT 


Mox+t 

x+M0t 


/Mq£-1  U 
\2 


+  arc  sin 


(18b) 


Region  C  (between  lines  x  =  t,  t  =  0,  and  x  =  c0-Mot) 

Ap  _  4a 
q  "  Mo 


(18c) 


From  the  pressure  distributions  it  is  possible  to  calculate  the 
indicial  lift  coefficient  Cj^t)  as  a  function  of  Mq  and  t. 
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Since 


=  ~f 

coJ 


Ap 

qa. 


dx 


the  following  results  are  obtained: 


First  time  interval  0  < t  < 


1+Mq 


Ct  (t)  =  — 
ija  M 


(19a) 


Second  time  interval  -CP  <t<  c‘ 


1+Mq  Mo-1 


cLa(t) 


-)*  -a— 


1  /n  c0— M  t' 

—  (  —  +  arc  sin  — - - i  +  —  arc  cos 

Mo  V2  t  /  VMo2  - 


t+M0c0-tM0; 


(19b) 


Third  time  interval  — — 2_<  t 
- -  Mq-1 


Ct  (t)  =  Z 


(19c) 


These  results  have  been  discussed  in  reference  7  for  values  of 
M0  greater  than  one.  They  still  hold,  however,  for  sonic  flight 
speeds  and,  in  fact,  can  be  reduced  to  the  expressions: 

First  time  interval  0  <t< 

-  2 
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Clte(t) 


=  4 


(20a) 


Second,  time  interval 


£<t 

2 


ciu(t)  =  «(l +  arc  sln  ^ +  <20b) 


The  indie ial  lift  coefficient  is  seen  to  be  constant  and  equal  to  4a 
up  to  the  time  t‘  =  —  or  up  to  the  time  required  to  travel  one-half 

2ac 

chord.  Following  this  first  time  interval  the  indicial  function  rises 
monotonically,  reaching  an  infinitely  large  value  as  time  increases. 

The  growth  of  Ct  (t)  is,  of  course,  in  agreement  with  the  fact  that 

the  steady-state  load  coefficient  becomes  infinitely  large  in  linear 
theory  for  M^l.  This  means  that  the  theory  cannot  be  used  to  predict 
the  complete  extent  of  the  Cj^(t)  variation  with  time  but  that 

during  the  earlier  part  of  the  motion  the  assumptions  remain  valid. 

In  figure  3,  curves  of 


<yt) 


2MJ 

c. 


are  plotted  as  functions  of 


for  values  of  MQ=1  as  given  by 
equations  (20a),  (20b),  and 
M^l.2,  Mq=1.4  as  given  by  equations 

(19a),  (19b),  and  (19c).  Also 
included  in  the  figure  are  variations 
of  Ct  (t)  for  Mo=0  as  calculated 
"^a 

from  Wagner's  results  (reference  8) 
by  E.  T.  Jones  (reference  9)  and 
also  for  Mg=0.8.  The  derivation 
of  results  leading  to  the  Mq=0.8 
curve  will  be  given  subsequently 
in  this  paper.  The  value  of 


Figure  3.—  Indicial— lift-curve 
slope  for  Mach  numbers  between 
0  and  1.4  shown  to  time  required 
to  travel  12  half -chord  lengths. 


0  (t) 


at  Mq=0.4  for  a  short  interval 
of  time  is  also  drawn.  The  dashed 
portions  of  the  curves  were  not 
calculated  but  were  drawn  to  agree 
with  the  known  asymptotic  value 
of  the  lift  function. 
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Application  of  the  indicial  lift  function  at  1^=1.—  Once  the 

indicial  lift  function  is  known,  it  is  possible  to  determine  the 
lift  corresponding  to  a  given  variable  motion.  Consider,  as  an 
idealized  example,  the  case 
where  the  airfoil  experiences 
an  abrupt  rising  and  sinking 
motion  at  regular  time  inter¬ 
vals  .  Such  a  motion  involves 
abrupt  plus  and  minus  angles 
of  attack  without  rotation  or 
pitching  so  that  a(t’)  is 
given  by  the  meander  or  square- 
wave  function  shown  in  figure 

4(a).  In  this  example  the  (a)  Impressed  angle  of  attack, 

variation  of  a  is  such  that 
the  curve  for  Cb( t  * )  can  be 

calculated  easily.  In  figure 
4(b)  Cj^t’)  is  shown  for  the 
case  in  which  the  discontin¬ 
uities  occur  at  intervals  of 
time  equal  to  cQ/V0,  that  is, 
after  each  chord  length  of 

travel.  The  principal  point  (b)  Eesulting  variation  of  lift, 

of  interest  in  this  example 
is  the  fact  that  such  a 

motion  yields  no  excessive  Figure  4.-  Lift  resulting  from 

value  of  lift  or  perturbation  square^vave  angle-of-attack 

velocities  and  the  entire  variation, 

analysis  is  within  the  frame¬ 
work  of  linear  methods. 

When  the  variable  motion  is  more  complex  in  character  the  lift 
coefficient  can  be  expressed  by  means  of  DuhamelTs  integral.  Corre¬ 
sponding  to  the  angle-of-attack  variation  a(t*)  as  a  function  of 
time,  the  lift  coefficient  Cj/t*)  is  given  by  the  expression 


j  p  t  * 

ClU’)  =  J  Citt(t,-T»)  a  (t«)  dT» 


In  analysis  related  to  equation  (21)  it  is  convenient  to  employ 
techniques  associated  with  the  use  of  the  Laplace  transformation.  (See 
reference  10.)  Thus,  if  the  Laplace  transform  f(s)  of  the  function 
f(t)  is  defined  by  the  relation 
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00 

f(s)  =  /  e"8*  f(t)  dt 


then  equation  (21)  can  be  rewritten  in  the  form 


Cl(s)  *  sCi^Cs )a(s) 


(22) 


Consider  now  the  case  of  a  lifting  flat  plate  oscillating 

harmonically  without  pitching  at  a  frequency  cs  and  maximum  angle 

of  attack  equal  to  a _  .  Setting 

max 


/.x  last  issut’ 

*(t)  =  ®  =  Onnx  8 


(23) 


then  equation  (22)  yields 


Pl(s)  _  4 

s— ico 


C/5 


co 


Jb 


-  erf .  /-*> 


(24) 


By  straightforward  manipulation,  the  inverse  transformation  of 
equation  (24)  can  be  shown  to  give 


-2l.  .  i*10*  (i  +  S 
aaax  l  V *v 


-iV 


C(cot— v)  —  i  S(o>t— v) 


-  2  i  Ni(cut,  v)  -  2y|  N2(o)t,  v)|  -  f  arc  tan/v/^^a  (25) 


where 
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co*  is  true  impressed  frequency 
io)lt f 

(«  =  “max  e  )• 

C(ayt— v),  S(cot— v)  are  Fresnel* s 
integrals  (reference  11). 

Ni(cot,v),  ^(cnt^v)  are  integrals 
defined  in  the  appendix. 

% 

If  the  response  to  a  cosine 
variation  of  a  is  required , 
only  the  real  part  of  equation 
(25)  is  used.  Such  a  response 
in  the  early  stages  of  the 
maneuver  is  shown  in  figure  5. 

For  very  large  values  of  time 
CL  as  given  hy  equation  (25) 

approaches  the  value 


Figure  5-—  Lift  resulting  from 
cosine-wave  angle— of-attack 
variation. 


CL  _  1  —  icut  (  f~i  iv.  , _ \ 

\,/i^ve  -erfTTvJ  (26) 

from  which  both  the  amplitude  and  phase  shift  of  CL  resulting  from 
either  an  impressed  sine  or  cosine  variation  of  a  can  be  readily 


Figure  6.-  Amplitude  of  oscillatory  lift  resulting  from  a  cosine 
angle-of-attack  oscillation  (without  pitching)  at  M  =  1. 


16 


MCA  TN  No.  1824 


determined.  Figure  6  shows  the  amplitude  of  lift  oscillation  corre¬ 
sponding  to  continuous  angle— of-attack  oscillation,  as  determined 

from  equation  (26),  plotted  as 


Figure  7.-  Variation  of  cos  cot 
with  t  for  various  values  of 
0). 


a  function  of  v 


It  is 


apparent  that  as  <x>  approaches 
zero,  CL  approaches  infinity; 


that  is,  as  the  impressed  wave 
approaches  the  "step"  function  ^ 
(fig.  7),  the  lift  coefficient 
approaches  infinity.  This 
result  is  in  agreement  with 
equation  (20h)  'as  t  approaches 
infinity.  As  the  frequency 
parameter  V  is  increased, 
however,  the  value  of  Ct  /% 


is  reduced  and  reaches  a  minimum  of  about  3-^  for  a  value  of  V  =0.9. 
For  a  speed  of  sound  around  1000  feet  per  second  and  a  wing  chord  of 
6  feet,  this  would  correspond  to  a  frequency  of  47.7  cycles  per  second, 
a  value  well  within  the  range  of  practical  flutter  frequencies.  It  is 
interesting  to  note  that  figure  6  also  shows  that  as  the  frequency  of 
oscillation  becomes  large  (i.e.‘,  v>3)  the  value  of  ^Lmax/ amax 

approaches  the  value  4.  This  is  the  same  as  the  value  for  C^t) 

in  the  early  stages  following  a  step  variation  of  a. 


Unsteady  State ,  MQ  < 1 

It  has  been  pointed  out  that,  in  the  determination  of  the 
indicial  lift  function  for  a  wing  traveling  at  subsonic  speeds,  the . 
lifting-surface  analogue  involves  the  calculation  of  load  distribution 
over  a  swept— forward  wing  with  subsonic  edges.  This  means  that 
recourse  cannot  be  made  in  the  solution  to  the  simple  source  distri¬ 
bution  method  used  in  reference  7  "to  treat  the  M<2=1  case.  Since, 
however,  a  portion  of  the  leading  edge  in  the  present  case  is  still 
supersonic,  the  problem  is  particularly, adapted  to  lifting-surface 
methods  developed  by  Eward  in  reference  12.  Figure  8  indicates,  as 
in  figure  1(b),  the  geometry  associated  with  the  boundary  conditions. 

In  the  Evvard  analysis  the  solutions,  as  in  the  previous  case,  are 
calculated  for  various  regions.  As  an  aid  in  identifying  the  different 
results  the  sketch  denotes  these  regions  by  Roman  numerals. 
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The  load  distributions 
and  indicial  lift  functions 
will  be  given  here  for  values 
of  t  up  to  the  time  when 
the  characteristic  from  the 
trailing  edge  first  crosses 
the  leading-edge  trace,  that 
is,  for 

1-Mq 


or 


<  Co 
a0(l-Mo) 


x 


This  period  covers  the  time 

in  which  the  wing  travels  Figure  8.—  Regions  used  in  the 

2M0/(l-MQ)  half  chords  and,  study  of  subsonic  unsteady  lift. 

since  the  present  analysis 

is  concerned  with  values  of  Mach  number  near  one,  will  in  some  cases 
extend  beyond  the  range  of  the  linear  theory. 


The  following  results  are  obtained  for  load  coefficient: 


Region  I  (between  lines  x=t,  t=0,  and  x=cQ— t) 

Ap  4a 

Region  II  (between  lines  x=-MQt ,  x=t,  and  x=cD-t) 


(27a) 


Ap  _  8a  /  Mp  /  t— x 
q  «Mq  Vl+Mo/  V+x 


Region  III  ^between  lines 


t  = 


x=cQ— t,  x=t,  x=— t  + 


g°o 

I+Mq’ 


(27b) 


and 
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U3  -4sl 
3  «% 


(v-v1)(u+v1) 


jsln-l  (m~Vl)(l-Mo)[(Y-l7i)(l+Mo)— C0iy/2  ]  +2vi [ v(  1+Mp )-Cq  V&-u(  1-Mq )  ]  1 
^  [v(1+Mq)— c0  J2  +  Vid-Mo)]  [uCI-Mq)  -  vi(l+MQ)]  J  1 

where 

U  =  ■—  (t-x) 

a/2 

▼  =  -t  (t+x) 

ys 

E(jJ’k)  }  incomplete  elliptic  integrals 
K,  E  complete  elliptic  integrals 

?(aJ  the  of  pressure  distribution  with  time  is 

shown  at  subsonic  speed  for  the  period  of  time  covered  by  equations  (27) 
For  purposes  of  comparison,  H  Kl ' 

pressure  changes  calculated 

from  equations  (18)  are  shown  ^  h  ^  t  /i  * 

in  figure  9(b)  for  supersonic  q 

flight  velocities.  ' 

Equations  (27)  suffice  \\  f  'dZj 

to  determine  the  initial  j  1  L — ^ 

growth  of  indicial  lift  coeffi-  !  j 

cient  at  subsonic  speeds.  j  raHpppft  \ 

Such  results  were  given  in  '  r  ”  ^jgjjjf 

figure  3  at  1^=0. 8  along  j0H0U 

with  the  calculated  growth  (a)  Subsonic.  (b)  Supersonic, 

for  about  one  chord  length 

rllZT  c  T  „  FlgUIe  9-  distribution  on 

value  of  CTi  (t  )  at  t  =  0  wines  receiving 


value  of  Cj^(t  )  at  t  =  0 

is,  for  all  flight  speeds, 
equal  to  4/M0* 


(a)  Subsonic.  (b)  Supersonic. 

Figure  9.-  Pressure  distribution  on 
wings  receiving  sudden  angle-of- 
attack  change  at  t  =  0. 


Expressions  for  C^(t )  are  as  follows 


First  time  interval  0  <t< 


1+Mo 
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Co-^U-«o) 


(28a) 


Second  time  Interval  -j-^-  <t  <  -^  ^ 


Cla(  '  '  nf^c0  l  1+MC  J  2t(l+M0)-ec„ 


+  It  \  t  “ 


(l+Mo): 


^t(l+Mo)-2c0]  [ 2c0— t ( 1-Mo2 ) ] 


k°*2ii ±5al  arc  sin/^fo  - 

1+MC  J  t(l+Mo)  1+Mo 


—  co“Mot 

+  1-0  -  arc  tanyg=}t  ±  f  - 

(28b) 


where  (m  is  given  "by  equation  (27e). 

tr 


PART  III  -  THREE-DIMENS IONAL  LINEAR  PROBLEMS  FOE  Mq  NEAR  ONE 


Steady  State 

General  solutions  for  arbitrary  Mach  numbers.-  Two  methods  of 

attack  are  available  for  the  solution  of  linearized  problems  at  sonic 
speeds.  In  the  first  place,  solutions  to  equation  (13)  can  be 
written  formally  and  the  extent  to  which  these  solutions  satisfy  the 
original  assumption  can  then  be  investigated.  In  the  second  place, 
general  solutions  of  equation  (7)  can  be  studied  in  the  limit  as  MQ 
approaches  1.  Since  this  latter  method  furnishes  added  information 
concerning  the  variation  of  the  variables  with  Mo,  it  will  be  used 
first . 
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In  linearized  theory  the  boundary— value  problems  of  wing  theory 
are  concerned  with  two  separate  properties  of  the  wing:  the  thickness 
effects  and  the  effects  produced  by  the  twist,  camber,  and  angle  of 
attack.  The  first  is  called  the  nonlifting  case  and  the  second  is 
the  lifting  case.  Solutions  of  equation  (7)  for  1  are  given 

in  reference  (13)  as  follows: 

In  the  nonlifting  case 


cp(x,y,z)  =  "Sc// 


Aw0(x1,yi)dx1dyi 
7  ( )  2-p2(  y-y  1 ) 2-02z2 


(29) 


where  (3  =A2-i  and  Aw0  =  2w0  where  wc  is  the  vertical  pertur¬ 
bation  velocity  on  the  wing  and  therefore  related  directly  to  the 
slope  of  the  wing  surface  relative  to  the  x  axis.  The  integration 
region  t  is  the  area  on  the  wing  within  the  Mach  forecone  from  the 
point  x,  y,  z. 


In  the  lifting  case 


cp(x,y,z) 


_1 
2  * 


p  p  02z  A9e(xi,yi)dxidyi 

4  J 


(30) 


where  AcpQ  is  the  jump  in  the  value  of  the  velocity  potential  in 
the  plane  of  the  wing.  The  sign  |  denotes  "finite  part"  of 

the  integral  and  introduces  special  integration  techniques.  (See 
reference  13.) 

Equation  (29)  expresses  the  velocity  potential  for  the  symmetri¬ 
cal  wing  in  terms  of  an  integral  involving  supersonic  source  distri¬ 
butions  while  equation  (30)  employs  doublet  distributions.  In  the 
two  cases  the  distributions  are  determined  from  the  geometry  and  the 
load  distribution  over  the  wing, respectively. 

Source  and  doublet  distribution  effectiveness  at  infinity.—  It 

is  well  known  that  the  lift,  drag,  and  pitching  moment  of  a  given  wing 
may  be  calculated  either  from  direct  integration  of  the  local  pressures 
on  the  wing  or  by  means  of  momentum  considerations  where  the  induced 
velocities  of  the  wing  are  determined  at  an  infinite  distance  and 
the  desired  forces  are  related  to  an  integration  over  a  control  surface 
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enclosing  the  wing.  In  the  three  following  sections  the  latter 
approach  will  he  considered  and  the  limiting  value  ,of  drag  at  MQ=1 
computed.  The  initial  portion  of  this  theory  requires  the  evaluation 
of  source  and  doublet  effectiveness  at  infinity  and  the  concept  of 
equivalent  source  position,  an  idea  which  appears  to  have  been  given 
first  by  W.  D.  Hayes  in  reference  l4. 

Consider,  as  in  figure  10,  a  point  P  with  coordinates  x,  y,  z 
lying  within  the  induced  field  of  a  supersonic  wing.  The  Mach  f ore- 
cone  from  P  is  given  by  the  relation 


x-xi  =  3v(y-yi)2  +  (z-zi) 


where  xi,  yi,  zi  are  running  coordinates  of  a  point  on  the  surface 

of  the  cone.  Introducing  polar 

coordinates 

y  =  r  cos  0,  z  =  r  sin  0 

and  rewriting  the  abscissa  of 
P  in  the  form 

x  =  x0+3r 

it  follows  that  the  trace  of 
the  forecone  in  the  zi=0 
plane  is,  in  the  limit  as  r 
approaches  infinity , 

xi  =  3yi  cos  0  +  xQ  (32) 

It  is,  moreover,  possible  to  show  that  the  effect  on  the  velocity 
potential  at  the  point  P  as  r  approaches  infinity  is  the  same 
for  all  points  (xx,  yx,  0)  for  which  xi-3yicos  0  =  constant.  The 
value  of  this  effect  is 


-  and  follows  from  the  asymptotic  evaluation  of  the  supersonic  source 
potential 


Figure  10 Coordinates  used  in 
study  of  supersonic  source. 
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£  ^8o)Tsuft„:-j^“  r  p  fvhe  ■«“  •* 

trace  to  (Xi-Py^cos  9  0  0)  tho  i  +  so^rce  shifted  along  the 

~is-  -  ~ 


cp  = 


2* V23r(Xo-Xlj 


(34) 


and  is  equivalent  to  the  potential  „+  v  »  X1_ 

The  induced  velocities  at  P  due  to  a  source  source  at  (xx,  0,  0). 
immediately,  for  arbitrary  Mn  Ld  for  T  (Xl*  Yl’  °}  follov 

gradients  of  <p  in  equations  (33)  and  (34)°  It^R0^6^  f^omthe 
that  equation  (33)  is  a  function  +,  11  is  ^Portant  to  note 

that,  in  general' ,  a  to^cTallT  f  t  *  azlmuthal  angle  of  P  so 
vith  respect  to  its  poS^iaHt  innMtV  Position 

is  independent  of  the  Jt££  ^  Ration  (34),  however, 

sy^i:^  16  to  the  study  of 

the  doublet  potential6  8  *  lifting  surfaces  are  to  be  analyzed, 


<P  = 


32z 

,2  ,2 


2*  [ (X-X1)2^|2(y_y1  )2_p2z2]3/2 


vith  rLpecrtrtL“otentiaieatiSf2-t<1UiV?lent  doublet  Position 
doublet  position  canavain  5*2?*?  ty  arises*  In  this  case  the 
Mach  cone  from  P  “  gSnSv  L  tt  to  the  trace  °f  the 
the  expression  7  aBd  the  potentlal  at  P  is  given  by 


cp  = 


32z 


2n[2pr( XQ-^i+^cos  e)]3/S 


(35) 


and,  for  Mach  numbers  near 


one 


cp  = 


2»t  [  20r  ( xQ— xx )  ] 3/2 


(36) 


Momentum  relations.-  The  vectorial  force  t  m  „ 

b°ay  ^  “  — 1  ^  the  ”  f™r™r 


NACA  TN  No.  1824 


24 


(37) 


where  vector  notation  is  used  and 

o  subscript  indicating  free-stream  condition 

p,p  local  static  pressure  and  density 
f  local  perturbation  velocity  vector 

For  the  purposes  of  the  /^Tlinearized  theory  and  the 

-sr  rsr:  w 

fac.  is  at  constant.  (See 

figure  11.) 


From  linearized  theory, 

=  1  - 
Po 


2  U 

*o  % 


and 


-  Pc  ■  -  Pofv  +  l<u2+v!!'f”a,]  +  5 


Figure  11.-  Surfaces  used  in  study 
of  momentum o 


The  end  faces  of  the  cylinder 

may  be  denoted,  as  in  the 

figure  by  I,  II,  the 
tigure,  uj  x,  7  rnhpn  in 

curved  surface  t>y  I  • 
supersonic  flow,  if  a  distri 
tion  of  sources  is  restricted 
to  a  region  downstream  of  7 
the  drag  D  on  the  ^dy  corre¬ 
sponding  to  the  source  ^istribu 
tion  is  given  by  the  expression 
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D  =  \  po  J  ( M02-1 )  U2+V*+W2 

where  vr  is  the  radial  component  of  the  perturbation  velocity.  Ho 
loss  in  generality  results,  moreover,  if  the  surface  II  is  moved 
infinitely  distant  downstream  and  the  radius  of  the  cylinder  is  made 
arbitrarily  large.  The  notation  II  and  III  will  henceforth  refer  to 
this  particular  conf iguration. 

If  the  drag  of  a  lifting  surface  is  to  be  calculated,  the 
surface  and  its  vorticity  wake  are  replaced  by  doublet  distributions 
and  in  that  case  the  integral  over  region  II  in  equation  (38)  is 
called  the  vortex  drag  of  the  body  while  region  III  yields  the 
wave  drag.  It  has  also  been  shown  (see,  for  example,  reference  17) 
that  the  vortex  drag  of  a  supersonic  wing  is  a  function  only  of  its 
span  load  distribution  and  is  equal  to  the  induced  drag  at  subsonic 
speeds  for  the  same  span  loading.  If  a  finite  nonlifting  body  Is 
considered,  each  of  the  velocity  components  in  region  II  is  attenuated 
in  such  a  manner  that  its  contribution  to  the  vortex  drag  is  zero. 

The  integration  over  region  III  again  provides  the  wave  drag  for  the 
nonlifting  body. 

The  combination  of  the  results  given  in  this  and  the  last 
section  provides  a  method  for  finding  the  wave  drag  of  an  arbitrary 
body.  The  first  step  is  the  determination  of  the  source— sink  or 
doublet  distribution  corresponding  to  the  body  and  then,  by  means  of 
the  principle  of  equivalent  positions ,  the  sources  or  doublets  are 
moved  to  the  x  axis.  The  wave  drag  is  then  calculated  from  equa¬ 
tion  (38)  once  the  induced  velocities  on  the  control  surface  are 
known.  In  the  next  section  the  wave  drag  will  be  written  in  a 
different  form  and  the  drag  at  sonic  speeds  will  be  investigated. 

This  analysis  will  also  provide 
some  insight  into  the  range  of 
validity  of  the  sonic  theory. 

Evaluation  of  wave  drag  as 
MU  approaches  one .—  In  order  to 
study  the  drag  of  a  symmetrical 
body  at  zero  angle  of  attack,  it 
is  convenient  to  consider  the 
general  expression  for  the  velocity 
potential  given  in  equation  (29). 
Introducing  first  the  transforma¬ 
tion  (fig.  12) 


Figure  12.—  System  of  axes  in 
transformation  equation  (39). 


^-Po  /  / 

j  j  in  j 


uvrrd@dx 


(38) 
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where 


6  =  -  yi  tan  \x  +  xi 
tj  =  n  sec  m- 


(39) 


tan  p  =  0  cos  9 


equation  (29)  becomes 


q>(xfy,z) 


1  P  P 

2lt  J  Jt  aA (x-Ki)2-0^(y-yx)2-02z2 


Since,  however,  it  has  been  shown  that  <p  evaluated  infinitely  far 
away  from  the  wing  does  not  change  if  a  source  is  moved  along  the 
line  ^constant,  it  follows  that  the  source  strengths  can  he 
integrated  along  these  lines*  The  second  integration  is  then  along 
Tj=o  where,  from  equation  (39)j  !=xi  and  the  value  of  the  potential 
at  an  infinite  distance  is 


—sf  f*"*™** 


Setting 


it  follows  that 


f(xi,p)  =  -  cos  p  J'  Aw0(xi,T))dn 


i  px~®T  f(xi,p)dxi 


/  ,  IP  Ax 

-  s’ /  Tjs 


*i)2-02r2 


(40) 


and  this  is  the  same  as  the  potential  for  a  body  of  revolution  with 
source  strength  per  unit  length  given  by  f(xi).  The  induced  velocities 
corresponding  to  the  potential  in  equation  (40)  are  found  to  be,  after 
first  integrating  by  parts  and  using  the  notation  cydxi  f(xi,p)  = 
f’(xi,p)  together  with  the  relation  f(0,p)  =  0, 
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and 


u 


Sep  =  _1 

3x  2n 


f  *Ui,n)dxi 
V(x-«i)2-p2r2 


vr 


dep 

¥ 


(x-^i)f  t(x1,|i)dx1 
'/(x-xx)2-(B2r2 


(41) 


(42) 


Asymptotic  values  of  the  velocity  components  for  large  values  of  r 
are  readily  seen  to  be,  after  first  setting  x=xD+Pr, 


u 


f ’(xi^Jdxx 


and 


f  ,(xx,|i)dxx 
VXq-Xi 


(43) 


(44) 


Equations  (4-3)  and  (44)  may  be  used  together  with  equation  (38) 
to  give  for  the  value  of  drag  the  expression 


f,(xx,q)dxi  p°  f*U2,n)dx2 
,/xq-xi  Jo  a/x0-x2 


(45) 


Assuming  that  the  body  is  of  finite  length  so  tfyit  f*(x)=0  for 
x>Z  reversal  of  the  order  of  integration  yields  the  relation 


\ 


/ 
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If  equation  (46)  had  "been  derived  for  a  body  of  revolution,  then 
f(x)  would  have  "been  independent  of  the  angle  M-  and  in  that  case 
the  expression  for  drag  would  reduce  to  the  form 


D  = 


J '  f»(xi)f,(x2)^|xi-x2|dxxdx2 


(46a) 


/  / 

This  expression  was  given  by  von  Karman  in  reference  15 • 

For  the  study  of  the  drag  of  a  lifting  surface,  consider  now  the 
general  expression  for  the  velocity  potential  given  by  equation  (30). 
The  doublet  distribution  occupies  in  this  case  both  the  wing  plan  form 
and  the  wake  since  the  jump  in  4*  exists  also  in  the  vortex  wake. 

By  use  of  the  transformations  in  equation  (39)?  equation  (30)  becomes 


cp  (x,y,z) 


3gz  P  r  cos  u  &p(6,T))d£dri 

2«  Jr  J  [(x-^x)2-32(y-yi)^-32Z233/2 


and,  exactly  as  in  the  case  of  the  source  distribution,  this  can  be 
reduced  to 


<P(x,y,z) 

Setting 


3fz 

2n 


I 


cos  (j.  dxi 


[(x-xi)2-P2y2-P2z2 ]s/2 


/ 


ACp(xi,ti)dT) 


I 

g(xi,|i)  =  cos  H  r  A<P(xx,'n)dq 


it  follows  that 


cp(x,r,0)  = 


32r  sin  6 


2k 


TX  Pr  g(xx,n)dxx 
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Integrating  by  parts  and  using  the  fact  that  g(0,|i)=0 


and 


u  = 


dcp  =  6111  0 
dx  2«r 


l 


X_Pr  P^gUxi^dX! 

[  ( x— xj. ) 2— (32r2  ]  3/e 


sin  0 
2«r 


r  ^  nx-(3r  (x_Xl)gi(Xl<ti)dXl 

lr  J0-  /(x^)2-?^2 


(47) 


* 


x-f3r  (x-xxjpgrg^xi^jdxi 
[  ( x— xi ) 2— P2r2 ] 3/2 


} 


(48) 


where  g*(xi,|j)  indicates  - —  g(xi,|j.).  Setting  x  =  x  +3r  and 

0x3.  0 

letting  r  approach  infinity,  the  asymptotic  expressions  for 
equations  (47)  and  (48)  become,  if  g’(0,|i)  =  0, 


(49) 

(50) 


The  relations  just  derived  may  "be  used  in  conjunction  with 

equation  (38)  to  give  the  wave  drag  of  a  lifting  surface.  This 
result  takes  the  form 


g*  *(xi,n)dxi 
J  X0-Xi 


g*  t(xg,|i)dxg 
yx0-x2 


(51) 


* 
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In  the  wake  of  a  lifting  wing  the  function  g’ ’(xi,^)  =  0  and  if, 
moreover. 


g*’(xi,|i)dxi  =  gf(c0,|i)  -  g*(0,n)  =  0 


reversal  of  integration  in  equation  (5l)  yields  the  simpler  expression 


It  is  possible  to  draw  some  general  conclusions  from  equations 
(46)  and  (52)  regarding  the  wave  drag  of  wings  and  bodies  of 
revolution  without  the  necessity  of  detailed  applications  to  particular 
configurations.  It  is  apparent  immediately  from  equation  (46a)  that 
the  wave  drag  of  a  body  of  revolution  at  zero  angle  of  attack  is 
independent  of  Mach  number.  This  conclusion  does  not  apply,  however, 
to  the  nonlifting  wing  since  the  distribution  function  f(x,[i)  in 
equation  (46)  contains  the  variable  p  which,  in  turn,  is  a  function 
of  both  6  ‘and  (3.  As  M0  approaches  one,  the  study  of  the  non¬ 
lifting  wing  is  divided  most  conveniently  into  two  parts,  depending  on 
the  behavior  of  f(x,p). 

Consider  first  the  more  general  situation  in  which  f(x,p)  is 
not  zero;  that  is,  the  case  in  which  the  number  of  sources  does  not 
equal  the  number  of  sinks  along  the  line  £  =  constant.  This  means, 
when  Mq  is  1,  that  an  unequal  number  of  sources  and  sinks  appear  in 
the  transverse  or  yz  plane  and,  if  equation  (46)  is  applied,  either 
a  finite  or  an  infinite  value  of  drag  can  result.  The  limiting  value 
of  drag  at  sonic  speed,  obtained  from  integrations  of  surface  pressures, 
was  given  by  Stewart  and  Puckett  in  reference  (l6)  for  several  wing 
plan  forms,  all  of  which  had  nonvanishing  values  of  f(x,p).  If  the 
pressure  distribution  is  calculated,  however,  the  local  pressure 
coefficients  are  seen  to  become  infinitely  large  as  sonic  speed  is 
reached,  even  for  the  body  of  revolution,  so  that  the  assumptions 
of  the  linear  theory  are  violated  and  the  reliability  of  the  drag 
predicted  by  equation  (46)  can  in  no  case  be  assessed  even  though 
the  predicted  values  remain  finite.  Equation  (43)  shows  also  that 
when  control— surface  methods  are  used  to  compute  drag  at  Mq  =  1, 
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the  x  conqponent  of  induced  velocity  increases  indefinitely  when 
f(x,fi)  is  not  zero  and  that  the  theory  is,  therefore,  no  longer 
consistent. 


In  the  very  special  second  case,  that  is,  when  f(x,|j.)  vanishes 
for  all  values  of  0,  the  analysis  just  presented  breaks  down  at 
equation  (40).  It  is  clear,  however,  that  in  this  case  there  are 
equal  numbers  of  sources  and  sinks  in  the  |  =  constant  plane  and 
the  behavior  of  the  flow  field  at  infinity  is,  therefore,  exactly 
the  same  as  that  which  would  have  been  produced  by  a  distribution 
of  doublets.  Equations  (49)  and  (50)  give  the  velocities  induced 
at  infinity  by  an  arbitrary  doublet  distribution.  These  induced 
velocity  components  are,  in  terms  of  (3,  one  degree  higher  than  the 
similar  components  for  the  nonlifting  case.  The  values  of  both 
u  and  vr  can  thus  be  expected  to  approach  zero  for  all  values  of 
Mq  as  r  approaches  infinity  for  any  flow  field  generated  entirely 
by  doublets  or  by  an  equal  number  of  sources  and  sinks.  It  follows 
then  that  the  linearized  theory  for  lifting  surfaces  (generated 
entirely  by  doublets)  and  for  bodies  with  thickness  distributions 
such  that  f(x,p)  vanishes  (generated  by  an  equal  number  of  sources 
and  sinks  in  all  |  =  constant  planes)  is  entirely  consistent  as 
M0  approaches  one  and,  in  particular,  for  Mq  equal  to  one.  This 
being  true,  it  follows  immediately  from  equation  (52)  that  the  wave 
drag  of  a  lifting  system  is  zero  at  sonic  speed. 


Thickness  solutions  at  Mq  =  1.-  A  swept-back  wing  of  constant 

chord  and  infinite  aspect  ratio  is  an  example  of  a  practical  aero¬ 
dynamic  shape  for  which  an  equal  number  of  sources  and  sinks  occur 
in  every  yz  plane.  (See  fig.  13.)  Consider  the  case  in  which  the 
wing  cross  section  is  diamond  shaped  with  a  slope  equal  to  \  in  a 
plane  normal  to  the  leading  edge.  Then,  in  a  transverse  plane, 
(section  BB  of  fig.  13)  wQ  equals  ±Vo^  cos  ,  the  minus  and  plus 
signs  applying,  respectively,  to  the  left  and  right  of  the  ridge  line. 
Accordingly,  the  solution  of  the  problem  can  be  written  in  terms 
of  a  distribution  of  sources,  thus 


2jt  J, 


x  cot 


cot  i! 


V0X  cos  t  In  t(y-yi)2+z2]  dyi 


«• 


x  cot  i|r 


cot  i| r 


V0\  cos  i  In  t(y-yi)2+z2]  dyx  (53) 
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The  value  of  8<P/8x  can  immediately  he  found  to  he 


from  vhich  it  is  apparent  that  as  r  =  v  y2+z2  becomes  infinitely 
large,  Scp/cfcc  approaches  zero.  In  the  plane  of  the  airfoil,  that 
is,  for  z  =  0,  8<p/dx  becomes 
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f  =  cos  » 
\  dx  JQ  it  tan  i|r 


(y  tan  : x—  ^ c0  sec  V)(y  tan  ty— x+  icc  sec  t)" 

(y  tan  ijr— x)2  -* 


and,  using  the  definition  for  pressure  coefficient,  Cp 
this  becomes 


2X  cos  t  , 

— 7 - ,  I* 

it  tan  >|r 


sin  t|r  —  x  cos  t|r 


1 


(55) 


Equation  (55)  can  be  derived  by  entirely  different  methods. 
Perhaps  the  most  direct  of  these  alternative  derivations  is  the  one 
introduced  by  R.  T.  Jones  in  reference  18.  The  general  statement  used 
in  that  report  is  that  the  component  of  translational  velocity  of  a 
cylindrical  body  in  the  direction  of  its  long  axis  has  no  effect  on 
the  motion  of  a  frictionless  fluid.  Hence,  the  pressures  over  the 
ving  shown  in  figure  13  are  the  same  as  those  over  a  wing  moving 
normal  to  a  free  stream  with  a  velocity  VQ  cos  \|r-  Using  the  Prandtl- 
Glauert  correction  to  the  thin  airfoil  solution  of  a  two-dimensional, 
diamond— shaped,  nonlifting  section  exposed  to  a  free  stream  with 
velocity  YQ  cos  i|r,  one  obtains,  for  Mq  cos  t<l. 


(56) 


where  wQ  is  the  vertical  induced  velocity  on  the  upper  side  of  the 
z  =  0  plane  and  x’  is  measured  normal  to  the  leading  edge.  If 
this  solution  is  referred  to  the  axial  system  of  figure  13  by  the 
trans  format ion 


x'  =  x  cos  i|r  —  y  sin  \|r 


and  the  integration  is  performed  after  taking  the  partial  derivative 
with  respect  to  x,  the  resultant  expression  for  pressure  coefficient 
is 


Cp - 


2\cos2  i|r 


Zn 


X- 


1  c 

2  ' 


cos  i|r— ysin  \|r 


M 


it  yi-M^cos2* 


(57) 
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At  sonic  speed  this  equation  reduces  immediately  to 


~  2\  cos  \I r 

°p  ■  -  rsr*  ,n 


cos  i|r  -  y  sin  t 


(58) 


which  is  identical  to  equation  (55)*  The  result  expressed  by  the  two 
equations  is,  of  course,  not  new.  The  significant  point  is  that  the 
same  variation  in  pressure  coefficient  was  obtained  by  two  widely 
different  avenues  of  approach  and  that  the  result  obtained  from  the 
particular  methods  applicable  to  sonic  speed  theory  is  in  agreement 
with  that  derived  from  more  conventional  analysis. 

Lifting-surface  solutions  at  Mp  =  1.—  It  should  be  mentioned  at 

this  point  that  Eobinson  and  Young  (reference  19 )  have  shown  by  means 
of  linearized  theory  that  supersonic  triangular  wings  and  subsonic 
elliptical  wings  of  the  same  aspect  ratio  have  values  of  lift— curve 
slope  which  approach  a  common  and  finite  limit  as  Mq  =  1.  The 
present  section  of  this  report  is  concerned  only  with  the  study  of 
lifting  surfaces  at  a  fixed  sonic  velocity  but  the  results  to  be 
obtained  are  in  agreement  with  the  limiting  values  of  reference  19- 


A  further  application  of  the  results  in  this  section  can  be 
made  to  the  case  of  very  low  aspect  ratio  wings  at  arbitrary  Mach 
numbers.  This  viewpoint  of  the  theory  was  first  presented  by 
E.  T.  Jones  in  reference  20  and  applied  to  triangular  wings  while 
in  reference  21  extension  was  made  to  include  pointed  wings  on 
slender  bodies  of  revolution.  This  duality  of  interpretation,  that 
is,  to  all  aspect  ratios  at  sonic  speed  or  low  aspect  ratios  at  all 
Mach  numbers,  applies  to  all  solutions  of  three-dimensional  problems 
obtained  from  equation  (13)*  In  the  subsequent  analysis,  attention 
will  be  confined  to  swept— back  plan  forms  of  lifting  surfaces  with 
pointed  vertices  and  thus  doublets  will  be  used  exclusively. 

In  application,  the  two  types  of  boundary  conditions  to  be 
considered  are  as  follows: 

1.  Boundary— value  problem  of  the  first  kind,  loading  specified.— 

It  is  given  that  AuQ  =  Uy— Uj  =  0  over  the  xy  plane  except  for 
the  region  occupied  by  the  wing  where  2uu  =  -2u^  =  AuQ  =  f(x,y), 
the  function  being  determined  by  the  specified  loading.  Over  all 
of  the  xy  plane,  the  imposed  conditions  are  AwQ  =  0. 
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2.  Boundary- value  problem  of  the  second  kind ,  surface  specified. 

Over  the  xy  plane,  the  imposed  conditions  are  Aw Q  =  0  every¬ 
where  and,  except  for  the  region  occupied  by  the  wing,  Au0  =  0. 

Over  the  region  occupied  by  the  wing  wG  =  wu  =  w 2  =  f(x,y)  where 
f(x,y)  is  determined  by  known  camber,  twist,  and  angle  of  incidence 
(The  delta  notation  again  indicates  the  jump  in  the  value  of  the 
variable  at  the  z  =  0  plane. 

Subscripts  u  and  l  indi—  y 

cate  conditions  on  the  upper  y=b(x)  A\ 

and  lower  surface,  respec—  1  y  =  bz(x) 

tively,  of  this  plane.)  / 


The  nature  of  the  differential 
equation  shows  that  the  value  of  9 
is  a  consequence  of  boundary  condi¬ 
tions  along  lateral  strips.  If,  as 
in  figure  14,  the  two  leading 
edges  are  given  by  the  expressions 
y  =  bi(x)  and  y  =  bs(x),  the 
velocity  potential  is  express¬ 
ible  in  the  form 


Figure  l4„-  Swept-back  plan 
form  with  curved  trailing 
edge . 


q>(x,y,z)  =  ■£- 


z  rb2  Aq>0(x,yi)dyi 


(y-yi)2+z 


If  the  boundary— value  problem  is  one  of  the  first  kind,  the 
general  expression  for  cp  follows  from  a  direct  integration  after 
noting  that 

nJ 

A9  (x,y)  =  /  Au  (x,yi)dyi  (60) 

0  Jbi 

Since,  moreover,  load  coefficient  Ap/q  is  related  to  Au0  by  means 
of  the  equation 


it  follows  that  the  velocity  potential  cp  can  be  found  for  any 
prescribed  load  distribution  of  a  given  plan  form.  The  value  of 
vertical  induced  velocity,  evaluated  at  z  =  0,  then  suffices  to 
calculate  the  twist  and  angle  of  attack  of  the  wing. 
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If  the  boundary-value  problem  is  one  of  the  second  kind,  the 
vertical  induced  velocity  is  given  on  the  wing  and  the  load  distri¬ 
bution  is  to  be  found.  In  this  case  the  use  of  equation  (59)  leads 
to  the  consideration  of  an  integral  equation.  Since,  however,  this 
integral  equation  is  a  common  one  in  aerodynamic  theory,  certain 
established  methods  may  be  applied  to  it. 

After  noting  that  A<pQ(x,y)  =  0  at  the  leading  edge,  integra¬ 
tion  by  parts  and  introduction  of  the  relation 


yields  for  perturbation  potential  the  expression 


<P  = 


AvQ(x,yi)arc  tan 


y-yi 

z 


dyi 


(6l) 


In  the  limit  as  z  approaches  zero  the  derivative  of  cp  with  respect 
to  z  reduces  to  the  form 


A  Avn(x>yj-)dyi 
2rt  J_b  y-7i 


(62) 


For  a  given  distribution  of  wQ  over  the  plan  form  of  the  wing, 
equation  (62)  represents  an  integral  equation  to  be  solved  for 
Av  (x,y)  subject  only  to  the  condition  that  the  Kutta-Joukowski 
condition  is  satisfied  at  all  subsonic  trailing  edges.  Once  AvQ(x,y) 
is  determined  it  follows  that 


(x,yx)dy! 

and 

Ap  =  _2_  6A9q 
q  Vc  8x 


(63) 

(64) 


In  the  present  report  the  solution  to  the  wing  plan  form  shown 
in  figure  15  will  be  presented.  The  value  of  A^q  which  satisfies 
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equation  (62)  is,  in  region  1 


A<P  =  — 2w  Jx2  tan2  0-y2 


(65) 


and  in  region  2 

Acp  =  — 2wc  x  tan  0 


E(lr0,k0)  -kn»2  F(*0,kJ 


o ,  o 


where  E  and  F  are  defined  in  the  appendix  and  where 


(66) 


,i,  .  /  x2  tan2  0— y2 

+  =  arc  s in/— - - - ^ 

j/  x2  tan2  0-ai2 


(67) 


kn*  = 


ai 


x  tan  6 


=  Vl-ko: 


(68) 


The  equation  y  =  ai(x)  of 
the  trailing  edge  for  which 
equation  (66)  is  valid  is 
given  by  the  formula 


ai  = 


V 


E0~V 


(69) 


which  expresses  ai  expli¬ 
citly  as  a  function  of 


\x  tan  0/ 


This  par tic— 


Figure  15.—  Dimensions  and  regions 
used  in  discussion  of  swept-back 
wings . 


x  tan  6t 

ular  choice  of  trail ing- 
edge  shape  was  used  to 
simplify  the  analysis.  The 
resulting  plan  form  approaches  a  constant-chord  wing  as  the  span 
increases.  The  variation  of  ai  with  x  is  given  in  figure  16;  and 
figure  17  shows  the  relation  between  aspect  ratio  and  span. 


The  loading  coefficient  is  given  in  the  two  regions  (defined  in 
fig.  15)  as  follows: 
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Figure  16.—  Graph,  showing  trailing— 
edge  position  of  the  swept— back 
wings  studied. 

Region  1 


Figure  17.—  Relation  between 
aspect  ratio  and  wing  semi- 
span. 


qa 


4x  tan2  6 
x2  tan2  0— y2 


(70) 


Region  2 


|=^tan0  [E(t0,k0)  + 

This  load  distribution  is  shown  in  figure  18  for  a  triangular  and  a 
swept-back  wing.  It  is  seen  that  the  loading  at  sonic  speed  bears 
a  close  resemblance  to  those  found  at  higher  Mach  numbers.  Two 
similarities  of  note  are  the  discontinuity  in  the  pressure  gradient 
at  the  Mach  wave  originating  from  the  trailing  edge  of  the  root 
chord  and  the  satisfying  of  the  Kutta  condition  only  where  the 
trailing  edge  is  subsonic.  The  lift  and  induced  drag  coefficients 
are  given,  respectively,  by 


x  tan  6 J  x2  tan2  9— y2 
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Cl  _  1  A  (  l of. 
a  tan  9  2  tan  0  \  sQ2 


Figure  18.  Pressure  distributions  for  triangular  and  swept— 
back  wings  at  M0  =  1. 


and 


CDj 

a2  tan  8 


A 

tan  9 


4 


(s0/c0  tan  8) 


(73) 


where 


These  coefficients  are  plotted  as  a  function  of  A/tan  9  in  figure  19. 
It  is  shown  that  the  values  of  Cp^/a.2  tan  9  and  Cj^/ a  tan  9  for 

finite  aspect  ratio  swept-back  wings  are  always  less  than  the  corre¬ 
sponding  values  for  the  triangular  wing  (A/tan  0=4). 


IL 


ko 
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When  the  span  of  the  swept—’ back  wing  becomes  very  large,  the 


A- » 

tanfl 


Figure  19. —  Variation  of  lift 

and  drag  with  aspect  ratio  for 
a  swept— back  wing  at  MQ  =  1. 


slope  of  the  trailing  edge 
approaches  asymptotically  the 
slope  of  the  leading  edge.  It 
follows  that  for  infinitely 
large  aspect  ratio  the  limiting 
value  of  the  load  distribution 
on  the  outboard  sections  should 
approach  the  value  given  by 
simple  sweep  theory  for  an 
infinitely  long  swept— back 
lifting  surface  with  constant 
chord .  This  result  is ,  in 
fact,  a  consequence  of  equa— 
tion  (71). 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
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APPENDIX 

LIST  OF  IMPORTANT  SYMBOLS 


y  coordinate  of  trailing  edge,  y  =  ai(x) 
free— stream  speed  of  sound 
local  speed  of  sound 
critical  speed  of  sound 


aspect  ratio 


( span) 


L  (wing  area) 
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free— stream  velocity 
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vr  radial  component  of  perturbation  velocity 

x,y, z  Cartesian  coordinates 

a  angle  of  attack  in  radians 

P  v/|i-m02| 

y  ratio  of  specific  heats,  for  air  7  =  1.4 

A90,Au0,Aw0  discontinuity  in  component  in  z  =  0  plane 

G  semivertex  angle  of  swept  back  wing 


V 


q>cp 

2 


p0  free-^tream  density 

$  total  velocity  potential 

cp  perturbation  velocity  potential 

o>f  impressed  frequency  (reference  to  true  time) 


0) 
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